Zero temperature transition from d-wave superconductor to underdoped regime 
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By using mutual flux-attaching singular gauge transformations, we derive an effective action de- 
scribing the zero temperature quantum phase transition from d-wave superconductor to underdoped 
regime. In this effective action, quantum fluctuation generated vortices couple to quasi-particles by 
a mutual statistical interaction with statistical angle 8, the vortices are also interacting with each 
other by long range interactions due to charge fluctuation. Neglecting the charge fluctuation, we 
find a fixed line characterized by 8 and calculate the universal spinon, vortex and mutual Hall drag 
conductivities which continuously depend on 8. Implications for double layer quantum Hall systems 
are given. When incorporating the charge fluctuation, we emphasize the importance of keeping the 
periodicity of the mutual CS term at 8 = ±1/2 to get the correct critical behaviors and point out 
possible future directions. The connection with Z2 gauge theory is discussed. 



I. INTRODUCTION 

In this paper, we are trying to study the nature of 
zero temperature quantum phase transition from d-wave 
superconductor at x > x c to the underdoped regime at 
x < x c of the high temperature superconductors ( Fig.l). 

T 
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underdoped regime d-wave superconductor 

Fig 1: The temperature ( T ) versus doping ( x ) diagram of 
high T c cuprates. 

It is well known that inside the superconductor phase 
x > x c , there are low energy quasi-particle excitations 
near the four nodes of the Fermi surface, the quantum 
phase fluctuation of the Cooper pair condensates are sup- 
pressed, the positive and negative vortices are bound to- 
gether. As the doping decreases from the right to x c , the 
phase fluctuation increases. At the underdoped regime 
x < x c , the quantum phase fluctuations are so strong 
that they generate free hc/2e vortices, destroy the long 
range phase coherence of the d-wave superconductor, 
therefore the superconducting ground state. However, 
the local short-range pairing still exists, the low energy 
quasi-particles at the four nodes remain. This quantum 
phase transition is driven by the condensation of hc/2e 
vortices. Near the transition, there are strong coupling 
between the quasi-particles and the phase windings of 



hc/2e vortices. 

A lot of work has been done on a closely related prob- 
lem where quasi-particles are coupled to the vortices gen- 
erated by external magnetic field inside the supercon- 
ducting state |(|. In the following, we will first review 
these work, then will extend the methods developed in 
these work to study this quantum phase transition where 
the quasi-particles are coupled to the quantum fluctua- 
tion generated vortices. 

Employing semi-classical approximation, Volovik 
pointed out that the circulating supercurrents around 
vortices induce Doppler energy shift to the quasi-particle 
spectrum, which leads to a finite density of states at the 
nodes Q. This effect ( Volovik effect) has been employed 
to attempt to explain the experimental observations Q| 
of longitudinal thermal conductivity K wx in . 

Starting from BCS Hamiltonian, Anderson || em- 
ployed the first single-valued singular gauge transforma- 
tion to study quasi-particle dynamics in the mixed state. 
Unfortunately, Anderson made an incorrect mean field 
approximation which violates the " Time-Reversal" sym- 
metry |(| , therefore leads to the incorrect conclusion that 
there is Landau level quantization of energy levels of the 
quasi-particle. Franz and Tesanovic (FT) employed a 
different single-valued singular gauge transformation to 
map the quasi-particle in a square vortex lattice state to 
Dirac fermion moving in an effective periodic scalar and 
vector potential with zero average and studied the quasi- 
particle spectrum numerically . They did not see the 
signature of Landau level quantization in their numerical 
calculations 

Unfortunately, in the vortex lattice state, the vector 
potential only provides a periodic potential instead of 
scattering quasi-particles, its effect on the energy spec- 
trum and other physical quantities is not evident. This 
motivated Ye to study the quasi-particle transport in the 
disordered vortex state where the random AB phase scat- 
tering may show its important effects. In Ref. M, Ye 
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observed that because infinite thin hc/2e vortices do not 
break the T reversal symmetry, any correct mean field 
theory should respect this T symmetry. He showed ex- 
actly there is no Landau level quantization due to this 
T symmetry. He applied the singular transformation in 
the FT gauge to disordered vortex state and found that 
the long-range logarithmic interaction between vortices 
suppresses the fluctuation of superfluid velocity (scalar 
potential), but does not affect the fluctuation of the in- 
ternal gauge field. Therefore the scalar field acquires a " 
mass " determined by the vortex density, but the gauge 
field remains " massless " . The quasi-particle scattering 
from the "massless" internal gauge field dominates over 
those from the well-known "massive" Volovik effect and 
the non-magnetic scattering at sufficient high magnetic 
field. This dominant scattering mechanism is a purely 
quantum mechanical effects which was completely missed 
by all the previous semi-classical treatments (2). In fact, 
it is responsible for the behaviors of both k xx and K xy in 
high magnetic field observed in the experiments . 

When the vortices are generated by quantum fluctu- 
ations or thermal fluctuations, themselves are moving 
around, new physics may arise. From Newton's third 
law, intuitively, the vortices must also feel the counter- 
acting AB phase coming from the quasi-particles. There- 
fore there is a mutual AB phase scattering between vor- 
tices and quasi-particles. In this paper, we enlarge the 
one singular gauge transformation developed in ||,^|]6| 
to two mutual flux-attaching singular gauge transfor- 
mations. We then apply them to study quasi-particles 
coupled to moving vortices near the zero temperature 
quantum critical point. In canonical quantization, we 
perform two singular transformations which are dual to 
each other to quasi-particles and moving vortices respec- 
tively. Just like conventional singular gauge transforma- 
tion leads to conventional Chern-Simon (CS) term, the 
two mutual singular gauge transformations lead to mu- 
tual CS term. This elegant mutual CS term precisely de- 
scribe the mutual AB phase scattering between vortices 
and quasi-particles. Alternatively, in path- integral pre- 
sentation, the effective action describes the quasi-particle 
moving in both vector and scalar potentials due to the 
phase fluctuations of quantum generated vortices. By a 
duality transformation presented in Refs. ||^] to a vor- 
tex representation, the quantum fluctuation generated 
vortices couple to quasi-particles by a mutual CS term, 
the vortices are also interacting with each other by long 
range logarithmic interactions due to charge fluctuation 
described by a Maxwell term. By carefully considering 
the periodicity of this mutual CS term, it can be shown 
that the effective action is essentially equivalent to the 
dual vortex representation derived by Senthil and Fisher 
(SF) from the elegant Z2 gauge theory except the 
Volovik effect missing in the Z2 gauge theory is also ex- 
plicitly incorporated into the present approach. Neglect- 
ing the charge fluctuation first and replacing the Z2 mu- 



tual CS theory by the U(l) mutual CS theory, we find 
a fixed line characterized by the mutual statistical angle 
and calculate the universal spinon, vortex and mutual 
Hall drag conductivities which continuously depend on 
9. This transition can also be viewed as a simplest con- 
finement and deconfincmcnt transition. The properties 
of the confinement and deconfincmcnt phases on the two 
sides of the critical point are also discussed. We stressed 
explicitly that the periodicity of Z2 mutual CS gauge the- 
ory on the lattice is not preserved by the U(l) mutual CS 
theory in the continuum. When the U(l) charge fluctua- 
tion is taken into account, we tentatively still replace the 
Z2 mutual CS theory by the U(l) CS mutual theory and 
treat both U(l) fluctuations on the same footing. We 
find the U(l) charge fluctuation turns the fixed line into 
a fixed point where charge ±e and spin 1/2 Dirac-like 
quasi-particle and charge 2e Cooper pairs are asymptot- 
ically decoupled. Although the result is pathological, it 
does suggest the spinon and holon are confined into elec- 
trons and Cooper pairs when condensing hc/2e vortices, 
in contrast to the condensing of double strength hc/e 
vortices discussed in ]T^ , |l2] | and reviewed in appendix C. 
Although we are unable to solve the critical behaviors 
of a combined Z2 and U(l) theory in this paper, we do 
stress that Z2 periodicity should be treated correctly to 
understand the critical behaviors of a theory with both 
discrete Z 2 and continuous U(l) gauge fields. 

The paper is organized as the following. In the next 
section, we derive the effective action with both the mu- 
tual C-S interaction and the charge fluctuation, first in 
the path integral language, and then in canonical quanti- 
zation representation. In Sec. Ill, we concentrate on the 
effect of the mutual C-S interaction and calculate the 
critical exponents and universal conductivities for gen- 
eral statistical angle 6. Implications for double layers 
Quantum Hall System are given. Some details are rele- 
gated to appendix A. In appendix C, we concentrate on 
the effect of charge fluctuation which is the only gauge 
fluctuation in double strength hc/e vortices. In Sec IV, 
we discuss the combined effect of mutual C-S interaction 
and charge fluctuation, we also stress the importance to 
keep the periodicity of Z2 gauge field when considering a 
theory with both discrete Z2 and continuous U(l) gauge 
fields. In Sec. VI, we discuss the connection of our ap- 
proach to the current fashion Z2 gauge theory and point 
out the possible open problems. In appendix B, we ap- 
ply the method developed in the main text to quenched 
random vortex array. 

II. THE EFFECTIVE ACTION 

Following the notation in Refs. we define d-\ — 

c-[(x),di = cj |Q. Because there are equal number of 
positive and negative vortices, we divide the positive vor- 
tices into two subsets 4> p i , <p P 2 and negative vortices into 
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two subsets 4>ni,4>n2- We define 4>a — 4>pi + <t>ni,4>B = 
4>p2 + 4>n2 and introduce the spinon by performing a gen- 



eral singular unitary transformation d — Ud s : 



U = 



a itpA 







,—i<j>B 



(1) 



where 4>a + <Pb = <j> is the phase of the Cooper pair. The 
original transformation is devised for static vortices. 
Here we extend it to moving vortices whose phase (f> is 
also fluctuating, therefore, depends on both the space and 
time. The spinon is charge neutral. In Anderson's gauge, 
4>a = 0, 4>b — i> or vice versa. In the former (latter), the 
spinon is electron-like ( hole- like). 

Expanding H s around the node 1 where p = (pf,0), 
we obtain the linearized quasi-particle Lagrangian C qp in 
the presence of the external gauge potential A^: 



r) + Vf(p x - a x )r 3 + v 2 (p y - a y )T 1 ]f/'i 



+ i>li>!VfV x (?) +ii>lT 3 <ipiv T (7 f ) + {1^2,x^y) (2) 



where v u 



\Afj_ is the gauge- invariant super- 



fluid momentum, it acts as a scalar scattering potential, 



bs) is the AB gauge field due to the 



phase winding of vortices |14|. Note that the external 
gauge potential only appear explicitly in the superfluid 
momentum v^. Because there are equal number of pos- 
itive and negative vortices in both <pA and 4>B> on the 
average, the vanishing of and is automatically en- 
sured, in this case, the Anderson gauge turns out to be 
more convenient than FT gauge. In the following, we use 
the electron- like Anderson gauge where a M = \d^§ [ p7| . 

We get the corresponding expression at node T and 2 
by changing Vf — > — Vf, v 2 — > — v 2 in the above Eq. 

£ l qp = ^\[{dr - ia T ) - Vf{p x - a x )r 3 - v 2 (p v - a a )r 1 ]^i 
- tpli/jiVfV x (r) + iiplT 3 ipiv T (r) + (1 2, x -> y) (3) 

Performing a P-H transformation ?/?i a = Ca^ip^ and 
the corresponding expression at node 2, the above Eq. 
becomes: 

£ l qp = ^i[( 9 r + ia T ) + v f (p x + a x )r 3 + v 2 (p y + a a )r 1 ]^i 
+ $l'ij)iVfV x {f f ) + i^^Vrir) + {1^2,x^y) (4) 

In order to make the final expressions explicitly SU (2) 
invariant, we perform the singular gauge transformation 
"012a = e~ % ^ A ~^ B ^^)\ a and the corresponding transfor- 
mation at node 2, then — > — a^, Eq.|] takes the same 
form as Eq.|| p5| . Adding the two equations leads to: 

£«P = WiaK^r - ia T ) + Vf(p x - a x )r 3 + v 2 (p y - a y )r 1 ]V'ia 
+ ip\ a ipiaVfVx(r) + itp\ a T 3 ^ la v T (f) + (1 -> 2,x -» y) 

where a = 1,2 is the spin indices. r' s matrices are act- 
ing on particle-hole space. As intended, the above Eq. 



is explicitly SU(2) spin invariant. Equivalently, we can 
start with the explicitly spin SU(2) invariant approach 
advocated in Ref. and perform the singular gauge 
transformation in the p-h space. 

C qp enjoys gauge symmetry U(l) u x U s (l) ( in 
fact, it is U(l) x Z 2 ), the first being uni- 
form and second being staggered gauge symmetry: 
Uniform (or external) U u (l) gauge symmetry 



<pA + X, <\>B -> <pB + X 



(6) 



Under this uniform U(l) transformation, the correspond- 
ing fields transform as: 

-> 4> + 2x, A a ^ A a + d aX 
V a -> v a , a a — > a a (7) 

d,v a ,a a all are invariant under this external U(l) 
transformation. Therefore the spinon d a is charge neutral 
to the external magnetic field. 

Staggered (or internal) U s (l) gauge symmetry 



4>a + x, 4>b -> 4>b-x 



(8) 



Under this internal U{1) transformation, the correspond- 
ing fields transform as: 



A a > A a 
a a -> a a + d a x 



(9) 



Although the spinon d a is charge neutral to the ex- 
ternal magnetic field, it carries charge 1 to the internal 
gauge field a a . 

It is easy to realize that U u (l) acts only on the boson 
sector, since the spinon is charge neutral, U s (l) acts only 
on the fermion sector. In fact, U s (l) should be a discrete 
local Z 2 symmetry because up and down hc/2e vor- 
tices are equivalent and do not break T symmetry ||. 

The phase fluctuation is simply 2 + 1 dimensional X-Y 
model: 



K 



K 



^ph = ^-<=-(d^-2A^) 



(10) 



After absorbing the scalar potential scattering part ( 
Volovik term) into £ p h, we can write the total Lagrangian 
H Cqp ~f~ Cph as. 

£ = 01a K#r _ ia r) + Vf(p x ~ O x )t 3 + V 2 (p y - a y )r 1 ]V'la 



(5) 



+ (1^2,^y) + |(^0-^) 2 



(11) 



where A e Jf = 2A }1 — K 1 J p and the quasi-particle elec- 
tric current is : Jo — x/jjT 3 ipj, J x = Vfip\ipi, J y = Vft/j^^- 
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From Eq.[ll], it is easy to identify the two conserved 
Noether currents: spinon current and electric current. 
The spinon current is given by: 



Jo = ipl{x)ipi{x) + ip\{x)ijj 2 {x) 

jy = iplfafa^ipiix) + ipl(x)v F T 3 i/j 2 (x) 



(12) 



Obviously the spinon current only comes from quasi- 
particle. In principle, the spinon current is not conserved 
due to scatterings between different nodes which lead 
to anomalous terms not included in Eqjll]. However, 
the inter-node scatterings involve large momenta transfer 
Ki — Kj for i j, so we neglect them due to momentum 
conservation in the long wave-length limit of the phase 
fluctuation |l(J . The exact conserved spin current jj* is 
with a/2 inserted in the above spinon currents. 

The electric current is given by: 



dC 
dA% ff 



= K(d u (f> - A e /f) = K{d a <j> - 2A^) + J u 

(13) 



Where the first part coming from Cooper pair and the 
second from the quasi-particle. Although they are not 
separately conserved, their sum is. 

The Noether current due to the the symmetry under 
4> — > 4> + X can be written as: 



3u 



dC 



= K(d u <f> A'Jf) ij' = 3l \jl (14) 



It is a combination of electric and spinon currents, 
therefore also conserved. 

Following Refs. we perform a duality transforma- 
tion to Eq.|l] 



K 



1 



Hid J + ifM - ijlA*Jf + ^-f eu 



(15) 



Where we have separated topological trivial spin-wave 
part and topological non-trivial vortex parts. 

Integrating out the spin- wave part, we get the conser- 
vation equation for the total current j* = j'Jj — ij*. In 
fact, as shown in the previous paragraphs, jfj and j" are 
separately conserved. Therefore we can introduce spin 
and electric gauge fields by: 



fn = e^xdual 
We can also define the vortex current: 



(16) 



Z7T 



where a" = d^(j)/2-K is the vortex gauge field. 

Substituting the above expressions into Eq.|lf 
reach: 

^/ e V + idrfjl - iA e /fe uvX d u at 



(17) 



we 



(18) 



where a„ 

h 1 



is the total gauge field felt by the 
moving vortices. 

Using $ for the vortex operator and adding the quasi- 
particle part, we get the following effective action: 

L = ipl[(d T - iaf) + v f (p x - a^)r 3 + v 2 (p y - a^ft^i 
+ (l-»2,a:-j/) 

+ \(d u - ial - ia^\ 2 + 7(|$|) + ^a% vX d u af 

+ \flv - iA e / f e uu \d u a\ - i", ,<),<!, (19) 

where 7(|$|) = m||<I>| 2 + g^l^l 4 + ••• stands for the 
short range interaction between the vortices. The last 
term is due to Berry phase in the boson representa- 
tion H which can be absorbed into A e Jf by redefining 
A e Jf — ► Afff + i^S^q, it acts like an external magnetic 
field in the z ( namely fi = ) direction. There are one 
species of vortex and N = 4 species of Dirac fermion( 
2 spin components at 2 nodes at the upper plane), the 
mutual statistical angle 9 = ±1/2. We also changed the 
notation by setting oft = a v , a* = a s ,a e = a. 

Although we derived the above equation in the 
electron-like Anderson gauge, in fact, it holds in any 
gauge. Eq.[H)] again enjoys the gauge symmetry U(l) x 
U(l), the first only acting on the vortex sector is the 
electric U(l) gauge field, the second on both the vortex 
and fermion sectors is the U(l) mutual CS gauge field. 
The mutual Chern-Simon term enforces the constraints: 
2 ^9j s u = e^ v xd v af,2n9j^ = e uv \d u a^. Physically, it 
means that when a quasi-particle encircles a vortex, it 
picks up a phase 2ir6. Equivalently when a vortex moves 
around a quasi-particle, it also picks up a phase 2n9. Al- 
though the conventional C-S term breaks T symmetry 
and has periodicity under 9^9 + 2 (l9|-|23) , the mutual 
C-S term does not break T symmetry and has periodicity 
under 9 —> 9 + 1. For example, 9 = —1/2 is equivalent 
to 9 = 1/2. In fact, this mutual CS term at 9 = 1/2 
should be the same as the Z 2 mutual CS term on the lat- 
tice discussed in Ref. |l(]] . Therefore strictly speaking oft 
and oft are Z 2 gauge fields, only a M is a U(l) gauge field. 
The Doppler-shifted term ( Volovik effect) is encoded in 
the last term in Eqjl^. This effect was not taken into 
account in the Z 2 gauge theory llO) . 
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In canonical quantization language, it is very instruc- 
tive to compare the well-known singular gauge transfor- 
mation leading to composite fermion in v — 1/2 system 
by Halperin, Lee and Read (HLR) |nj to the singular 
gauge transformations performed in the spinon-vortex 
system in this paper. The crucial difference is that in the 
former, we attach electron's own 9 flux to itself (9 = 2 
in v = 1/2 system to keep fermion statistics intact ), the 
constraint V x a(r) = 2n9p(r) in Coulomb gauge di<ii = 
leads to composite fermion i\) c coupled to conventional 
CS term which breaks T symmetry, has periodicity un- 
der 9 — > 9 + 2; however, in the latter, we attach vortex's 
9=1/2 flux to quasi-particle or vise versa, the two con- 
straints V x fff'(f) = 27r0j£(r),V x a*(r) = 2ir9j^(r) 
in Coulomb gauges diaf = diaf = naturally leads to 
quasi-particles and vortices are coupled by mutual CS 
term in additional to the electrical gauge field which me- 
diates the logarithmic interaction between vortices. The 
mutual CS term does not break T symmetry, has period- 
icity under 9 — > 9 + 1. 

It is important to point out that when a vortex is mov- 
ing around a closed loop, it pick up two phases, one is Z 2 
phase due to the spinon current, described by the mutual 
C-S term, another is U(l) phase due to the total electric 
charge current described by the Maxwell term. In the 
following section, we will neglect the electric charge fluc- 
tuation and concentrate the mutual statistics term. In 
Sec. IV, we will study the effect of charge fluctuation 
described by the Maxwell term. For completeness, we 
relegate the review of the effect of charge fluctuation to 
appendix B. 

For simplicity, we take the relativistic form for both 
fermion and boson, because the anisotropies in Eq.|l^ are 
expected to be irrelevant near the zero temperature QCP 




III. THE EFFECT OF MUTUAL STATISTICS 

In this section, we neglect the charge fluctuation, 
namely setting a M = in Eq.[H| The charge fluctuation 
can be suppressed by condensing hc/e vortices reviewed 
in appendix C. We could add the kinetic and potential 
terms for the hc/e vortex operator <E>2 to Eq.[l9|: 

U 2 = \(d^-i2a^)<S> 2 \ 2 + V(\<S> 2 \) (20) 

As shown in appendix C, there is no mutual statistical in- 
teraction between spinon and hc/e vortices $2- The long 
range logarithmic interaction between hc/2e vortices $ 
and hc/e vortices $2 is mediated by the electrical gauge 
field a M . Condensing < $ 2 >= $20 will generate a mass 

term -^-(a^)^ which dominates over the Maxwell term. 
Integrating out leads to 

^rl\(f»;") 2 - iA^e^M + (jl) 2 t \ (21) 



All the generated terms only renormalize the short 
range interactions already included in T^(|$|). The model 
is also closely related to the Hall drag problem in double 
layer Quantum Hall system ^6|, therefore the model it- 
self is interesting on its own right and deserves detailed 
investigation. 



A. Quantum Critical point 

In order to calculate conductivities, we add two source 
fields and A* for the quasi-particles and vortex re- 
spectively: 

+ \(d, - za* - zA*)$| 2 + V(m) + ^a%„ x d„af (22) 

where a = 1,2,3,4 stands for N = 4 species of Dirac 
fermion. In fact, two Maxwell terms for and a* can 
be added to the above equation, but they are subleading 
to the mutual CS term in the low energy limit. It is ex- 
pected that there is no periodicity under 9 — > 9 + 1 in 
the continuum limit. 

If the theory is properly regularized on the lattice, it 
should has the periodicity under 9 —> 9 + 1, for exam- 
ple, 9 = —1/2 should be equivalent to 9 = 1/2. In fact, 
0=1/2 mutual U(l) CS theory could be the same as the 
Z 2 mutual CS theory developed in p0[. 

The RG calculation in Refs. |20|-|22]| can be used to 
show that 9 is exactly marginal, therefore there is a line 
of fixed points determined by the mutual statistical angle 
9. In order to calculate the spin conductivity along this 
fixed line, a source field could be introduced to couple to 
the spin current = 4>t,Jt*(&)ai34>f3- Similar calculations 
follow. 

Integrating out both fermion and boson leads to: 

C = -\a%-k)n%(k)at(k) - \a${-k)TS${k)a*{k) 
- \at(-m%(k)at(k) ~ l -al(-k)Ult(k)ai(k) 
- Al(-k))e, vX k x (at(k) - At(k)) 
^«(- k ) ~ Al(-k))e^ vX k x (at(k) - At(k)) (23) 

where the exact forms of n' s are dictated by gauge in- 
variance and Furry's theorem : 

U^(k)=U 1 (k)k(S, v ~^-) 

n^(k) = u^(k)=n 2 (k)e^ x k x 

n%(k)=^(k)k(5^-^) (24) 



5 



Where IIi,n2,Il3 are the polarizations for 
a^a^,a t l'a' s ',a' s 'a' s '. 

If we are only interested in DC conductivities, for sim- 
plicity, we can put k = 0, Eqj23| becomes ( for the most 
general form, see appendix A): 

£ = --af(-uj n )Ilx\u) n \af (w„) - ^af (-w„)n 2 e i:) cj„a*(cj„) 
~ T^af {-uj n )n 3 \uj n \af (u n ) - ~af(-uj n )n 2 eijUJ n af(L> n ) 
- ^(af - (-cj n ))e i:j uj n (af(uj n ) - Af(w n )) 



o 



Fig 2: The renormalized propagators for eft (thick wiggle line ) 
and a* (thick dashed line). The thin wiggle line stands for c$ , the 
dashed thin line for a*, the thin solid line for fermion propagators, 
the thick solid line for the boson propagators. 



By using the bare propagators < 



4°* 



Ye {af{-u n ) - ^(-^eyc^o/K) - (25) 
If we define af(uj n ) = eyo*(w n ), Af (w n ) = eyA?(u/„), 



the above equation becomes diagonal in the spatial in- 
dices i = 1,2. Finally, integrating out af,af leads to 



L ~ 2 ( 1 ' 1 ' [a H a* ) { A* 



(26) 



Where spinon, vortex and mutual Hall drag conduc- 
tivities are: 



.1. 



Hi 



e' n 1 n 3 + (i/e-n 2 ) 5 



n. 



ITII3 + (i/6» - n 2 



a H = ( i )2 n2-0(nin 3 + ni) 



^ n 1 n 3 + (i/e-n 2 ) 2 



(27) 



In fact, all the three conductivities can be written in 
the elegant connection formula: 



Pij = (pfb) 



(28) 



~£fj,v\k\/k and the bare fermion and boson loop results 
n/," = 11^° = -f^fc(<^„ - k^k^/k 2 ), we can find easily 
the renormalized a* and a* propagators G^ v = G* y = 

— iQ^i^fiu ~ k^kv/k 2 ). In contrast to the conventional 
CS theory studied in 20[ 2l|] , the propagators are even in 
fc, this is because the theory respects T symmetry. On 
the other hand, in contrast to the Maxwell propagators, 
they behave as 1/fc instead of 1/k 2 . 

The three loop diagrams for II i are given by: 




WW 



Fig 3: The three loop diagrams of Eli. The thick wiggle line 
stands for the renormalized propagators of a*, the thin solid line 
stands for the fermion propagator. The one loop diagram is not 
shown. 

By using G^ v and extracting the symmetric part of the 
gauge propagator in the large N results in Refs. p0| , pl"| ], 
we are able to calculate the above three loop diagrams. 
Furry's theorem can be used to eliminate large number 
of null diagrams. We get the following series: 



with the conductivity tensor of fermion and boson given 
by 



o~fb 



if -n 2 
n 2 n 3 



(29) 



Although its form is similar to the conventional con- 
nection formulas discussed in [fl9|-pT||, the physical inter- 
pretations of the conductivities are quite different (see 
the following). 

When the vortices are generated by external magnetic 
field and pinned by impurities as discussed in jq] , the to- 
tal conductivity is the same as the fermion conductivity 
because the static vortex do not contribute. As explained 
in appendix B, the static vortices only feel the electric 
gauge field a M , but not the statistical gauge field a*. 

The above expressions are exact, but ni,n 2 ,II 3 can 
only be calculated perturbatively in the coupling con- 
stant g 2 — 2ir9. The renormalized propagators for oft 
and a* can be found from the following Feymann dia- 
grams: 



8 V 16 (2tt) 2 y 



(30) 



where N = 4 is due to the sum over 4 species of Dirac 
fermions. 

The three loop diagrams for II 3 are given by: 

o • o • o 



+ 



o 



Fig 4: The three loop diagrams of II3. The thick dashed line 
stands for the renormalized propagator of a*, the thick solid line 
stands for the boson propagator. The one loop diagram is not 
shown. 



For bosons, only one loop result is known: 
n 3 = Nl(l+g 4 + g 8 + ---) 



(31) 
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Note that both III and II3 are even functions of 9. 

From Furry's theorem, one of the first non-vanishing 
diagram for II2 is: 





Fig 5: One of the four loop diagrams of 1I2. the thin solid line 
stands for the fcrmion propagator, the thick solid line for the boson 
propagator. All the other four loop diagrams can be obtained by 
shuffling the positions of the three bare propagator lines. 



N(g 6 



which is an odd 



The series is IT2 
function of 9. 

From Eq.p7j, er^cr® are even functions of 9, but the 
mutual Hall drag conductivity is an odd function of 9. 
These are expected from P-H transformation. Under 
the P-H transformation of the vortex operator $ — > $t 
in Eq.E2, it can be shown that ^4* , 9 is equivalent to 



-A, 



From Eq.Eq, we reach the same conclusions. 



Specifically, a H takes opposite values for 6 = ±1/2, the 
periodicity under 9 — > 9 + 1 is not preserved. Exper- 
imentally, the Hall drag conductivity can be detected 
by measuring the transverse voltage drop ( or transverse 
temperature drop for thermal conductivity ) of spinons 
due to the longitudinal driving of vortices. The Hall 
drag conductivity in double layer Quantum Hall systems 
has been investigated by several authors pq . In double 
layer systems, the electrons in different layers are treated 
as two different species. There is a mutual CS interac- 
tion between the two species (both are fermions) which 
is directly responsible for this Hall drag conductivity, al- 
though the Coulomb interaction between the two species 
is responsible for the Coulomb ( longitudinal ) drag [^7j . 
This example shows that no external magnetic field is 
needed to produce a Hall effect ! However, 9 = ±1/2 are 
very special points, J7(l) mutual CS term reduces to Z2 
field, there should be no Hall drag conductivity just like 
a = ±1/2 vortex leads to no Hall conductivity. 

Again, by using the large N result of |2^,|2l|, we find 
the anomalous dimensions of the fermion and vortex to 
two loops: 



Fig 6: The two loop diagrams for the self energies of spinon (a) 
and vortex (b). The renormalized operators are used. 

The results are: 



48tH 



9 



12tt 2 



N 



(32) 



where t]xy ~ 0.038 is the anomalous dimension for the 
3d XY model ||. 

The correlation length exponent can also be calculated 
to two loops by the insertion of the operator <f>^<f>: 



Fig 7: The cross stands for the operator insertion of "J^c^. 



The result is: 



12tt 2 



N 



(33) 



where vxy ~ 0.672 is correlation length exponent for the 
3d XY model E3. 



It is instructive to go to dual representation of Eq,22 
namely go to the boson representation: 



1 



where V{\^\) = m 2 M\ 2 + + 



2 I JJ2 



d v {al 



/9)-iA% vX d v a{ (34) 



Integrating out a® leads to the constraint up to a pure 



gauge: 



,1> - 



Substituting the above constraint to Eq| 



(35) 
and setting 



dp, we find: 



£ = iPbn(dii - i0a n 
+ |(9 M -i<v)<£| 2 + 



f 2 



V(|,/>|) - iA%„ x d„a x (36) 



The above Eq. indicates that fermions and bosons are 
coupled to the same gauge field whose dynamics is de- 
scribed by Maxwell term instead of C-S term. This Eq. 
is simply 2 + 1 dimensional combination of spinor QED 
and scalar QED. The RG analysis at 4 — e by dimen- 
sional regularization is possible, because the marginal 
dimension of all the relevant couplings are 4. However 
a RG analysis directly at 2 + 1 dimension is formidable, 
the physical meaning and the exact marginality of 9 is 
obscure in the boson representation. However they are 
evident in the dual vortex representation Eq.|22|. It is 
evident that there is no periodicity under 9 — ► 9 + 1 in 

If we perform duality transformation again on Eq.36 
to go to the vortex representation, then we recover Eq. 19 
upon neglecting the two Maxwell terms relative to the 
mutual C-S term. 

In the next subsection, setting the two source terms 
vanishing, we look at the properties of the different 
phases on the two sides of this quantum critical point. 
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B. Disorder and superconducting phases 

In the disordered phase, the vortex condense < $ >= 
$o which generates a mass term for a* in Eq]22] 



<t> 2 



(37) 



where the subscript t means transverse projection. 



Integrating out the massive a* leads to a Maxwell term 



for at: 



(38) 



This is simply 2 + 1 dimensional spinor QED which 
was studied by large N expansion in |25| . In fact, we 
reach the same description from the boson representa- 
tion Eq.^6|. Because in the disordered phase, the boson 
<j> is massive, therefore can be integrated out, it generates 
the Maxwell term -r^— f, 2 „ which dominates over the ex- 
isting non-critical Maxwell term. We reach Eq.pS after 
identifying ~ *q. 

In the superconductor phase, the vortex $ is mas- 
sive, therefore can be integrated out, it leads to the old 
Maxwell term 4^(/*„) 2 - Integrating out a* generates 
a mass term: 



(39) 



where m$ is the mass of the vortex. 

The Dirac fermions become free. In fact, we reach the 
same description from the boson representation Eq.36 
Because in superconductor phase, the boson <f> condense, 

therefore generates a mass term 2 O^)* wmcn renders 
the Maxwell term ineffective. We reach the same conclu- 
sion from both sides by identifying m$ ~ 

In short, in the disordered phase, the system is de- 
scribed by spinor QED Eqj3^; in the superconductor side, 
by free Dirac fermion. We can view the transition as 
the simplest confinement-deconfinement transition. In 
the confined ( disordered ) phase, the boson and fermion 
are confined together by the fluctuating gauge field. In 
the deconfincd ( superconductor) phase, the boson con- 
densed, the fermion becomes free. There is a line of fixed 
point governed by the mutual statistical angle 9 separat- 
ing the two phases. 



term in the continuum. In this section, we try to con- 
sider the combined effects of Z 2 mutual statistical gauge 
fluctuation and U(l) electrical gauge fluctuation. After 
replacing the Z2 discrete gauge theory by a U(l) mutual 
CS theory, we treat 17(1) mutual CS theory and U(l) 
electrical gauge theory on the equal footing. Although 
the replacement leads to somewhat misleading conclu- 
sions, the heuristics derivation suggests that the conden- 
sation of hc/2e vortex condensation indeed leads to the 
confinement of spinon and chargon into Cooper pair and 
electron, in contrast to the condensation of hc/e vortex. 
We do not intend to understand the critical behaviors of 
a combined Z2 and U(l) theory in this paper. We simply 
stress that Z2 periodicity should be treated correctly to 
understand the critical behaviors of a theory with both 
discrete Z2 and continuous J7(l) gauge fields. 



A. Quantum critical point 



Putting a* — » a* — a M in Eq.[[{| we get 

C = ^(d, - iaf)^ + \(d,- ia*)*| 2 + 7(|*|) 

' <e^ A cU* - + A e / f ) + 1/^ (40 ) 



2tt9 m 

Integrating out the electric gauge field a M leads to: 

£ = ^7/*@m ~ ia+Wa + ~ ia*)*| a + V(l*l) 
1 



+ ^i^d v at + -{at/9 + Alfff t 



(41) 



Comparing to Eq.|2^, it is easy to see that the charge 
fluctuation leads to a mass term for the gauge field ajj. 
Shifting a^/9 + Afff — > a M and adding the gauge fixing 
term ^-(c^a^) 2 , we can integrate out the massive gauge 
field a M in Lorenz gauge a — and find: 

C = iph,(d, - < J )4 + - z<)$| 2 + 7(|*|) 
+ \(fp 2 - i(A e / f + Oj'^xduat + (iP? (42) 



Note the Maxwell term for a* is generated by the in- 
tegration over the massive o p . 

Setting Afj, = 0, integrating out the fermions only leads 
to higher derivative terms than the Maxwell term: 

C = \(d, - za*)$| 2 + y(|*|) + -(ft? + ■■■ (43) 



IV. THE EFFECT OF CHARGE FLUCTUATION 

In this section, we try to investigate the effect of charge 
fluctuation on the fixed line characterized by the statis- 
tical angle 9 discussed in the last section. In the last 
section, we pointed out the periodicity of Z2 gauge CS 
term on the lattice is not preserved by U(l) gauge CS 



where • • • means higher than second order derivatives. 
Therefore the vortex and fermion are asymptotically de- 
coupled. It indicates that the charge fluctuation ne- 
glected in the last section destroy the fixed line char- 
acterized by 9. However, very different conclusions are 
reached in Ref. pj| . We believe that the authors failed 
to treat the charge gauge field fluctuation correctly. 
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Just like the last section, it is instructive to go to dual 
representation of Eq.^l], namely go to the boson repre- 
sentation: 

C = - ia^ a + |(^ - la^ 2 + 



+ l Ui v ) 2 + iale flvX d v {a{-ai/e) 



(44) 



Integrating out a* leads to the same constraint as 
Eq|35| up to a pure gauge: 



at = Oaf, 



(45) 



Substituting the above constraint to Eq.|4| and setting 
af = a^, we find: 



£ = i'll^ - i9a^ a + - ia^\ 2 + Vm 



(46) 



Comparing to Eqj3^, the only difference is that the 
gauge field acquires a mass due to the charge fluctuation 
which renders the Maxwell term ineffective. Up to irrele- 
vant couplings, we can safely set a M = A e J* in the above 
equation and find 

£ = 4>ll»(d^ - i29A^ a + \{d, - i2A^\ 2 + V(\cf>\) + ■ ■ 

(47) 

where • ■ • means the irrelevant couplings between bosons 
and fermions. 

The above Eq. leads to in-consistent conclusion that 
the original charge neutral spinon tp becomes charged: 
carrying charge e for 9 — 1/2 and — e for 9 = —1/2. 
Therefore we also reach ambiguous conclusions for 9 = 
1/2 and 9 = —1/2. The reason is that we replace the Z 2 
CS gauge theory on the lattice which explicitly respects 
the equivalence between 9 = 1/2 and 9 = — 1/2 by the 
U(l) CS theory which does not respect this equivalence 
in the continuum limit. Identifying which is a discrete 
Z2 gauge field with which is a U(l) gauge field is a 
direct consequence of this problematic replacement. 



B. Disordered and superconducting phases 

We follow the discussions in the previous section. In 
the disordered phase, the vortex condense < <f> >= $0 
which generates a mass term for a* in Eq.^. Integrating 
out the massive a* leads to a Maxwell term for af.. 



^h^d^-iaf^ipa 



4<Dg 



M/9 + A<*f)* 
(48) 



The Dirac fermions become free and carry spin 1/2 and 
charge 29 e. In fact, we reach the same conclusion from 
the boson representation Eq.[f6|. 

In the superconductor phase, the vortex <!> is mas- 
sive, therefore can be integrated out, it leads to the old 
Maxwell term j^—{f^) 2 - Integrating out a* generates 
a mass term for af,: 

£ = th^d* - %4)1> a + ^(at) 2 t + \{a+/6 + Af/f) 2 

(49) 

where m$ is the mass of the vortex. 

Diagonizing the last two mass terms leads to a contin- 
uously changing charge . In fact, we reach the same de- 
scription from the boson representation Eq.^. Because 
in superconductor phase, the boson condense < <f> >= 'Po, 
therefore generates a mass term for a M : 

£ = ^7m(^ - iOa^a + f(a,) 2 t + \{a, + A**)* 

(50) 

Which is essentially the same as Eq.|^| after the identifi- 
cation m$ ~ (j>Q. 

Again, the pathological claims are due to the artifact 
of replacing the Z2 mutual statistical gauge theory by 
a (7(1) mutual CS theory, then treating Z2 gauge field 
and U(l) electrical gauge field on the same footing. Al- 
though all these results are the artifacts of this replace- 
ment, they do suggest that the mutual statistical interac- 
tion between hc/2e vortex and spinons leads to the con- 
finement of spinon and chargon into electron and Cooper 
pair. 



V. DISCUSSIONS AND CONCLUSIONS 

The original singular gauge transformation was 
proposed for static vortices. In this paper, we extend the 
singular gauge transformation to moving vortices gener- 
ated by quantum fluctuation. By making a close anal- 
ogy to the conventional singular gauge transformation of 
FQH system, we perform singular gauge transformations 
attaching flux of moving vortices to quasi-particles or vice 
versa. Just like conventional singular gauge transforma- 
tion leads to conventional CS term, the two mutual sin- 
gular gauge transformations lead to mutual CS term. In 
this way, we propose a intuitive and physical transparent 
approach to bring out explicitly the underlying physics 
associated with the condensation of the hc/2e vortices. 

Based on the earlier work by Balents, Fisher and Nayak 
p3[ , Senthil and Fisher developed Z 2 gauge theory |U0 11 
to study quasi-particles coupled to vortices generated by 
quantum fluctuations. By breaking electrons and Cooper 
pairs into smaller constitutes: chargons with charge e, 
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spin and spinons with charge 0,spin 1/2, SF intro- 
duced a local Z 2 gauge degree of freedom to constrain 
the Hilbert space to the original one. The effective ac- 
tion describes both chargons and spinons coupled to lo- 
cal fluctuating Z 2 gauge theory with a doping depen- 
dent Berry phase term. By the combination of standard 
duality transformation of 3 dimensional XY model and 
that of Z 2 gauge theory the action is mapped into a 
dual vortex representation where the hc/2e vortices and 
spinons are coupled by a mutual Z 2 CS gauge theory. As 
usual vortices in XY model, the hc/2e vortices also cou- 
ple to a fluctuating U(l) gauge field which mediates the 
long-range logarithmic interaction between the vortices. 
Starting from the dual representation, the authors in 
Rcf. |ll| studied a transition from d-wave superconduc- 
tor to confined Mott insulator driven by the condensing 
of hc/2e vortices. In order to study the critical behaviors 
of this particular confinement and de-confincmcnt tran- 
sition, they replaced the Z 2 mutual CS theory on the 
lattice by U(l) mutual CS theory in the continuum and 
performed renormalization Group (RG) analysis. Un- 
fortunately, some of their RG analysis are incorrect as 
demonstrated in Sec. IV. 

It was shown in the text that our resulting general ef- 
fective action Eq.|l9| at the two special statistical angles 
9 = ±1/2 is essentially equivalent to the dual vortex rep- 
resentation of Z 2 gauge theory developed by SF, except 
it also bring out explicitly the Volovik effect which was 
missed in SF's Z 2 gauge theory. 

As explicitly stressed in this paper, the £7(1) mutual 
CS theory in the continuum does not have the required 
periodicity under 6 — > 6 + 1. In order to preserve this pe- 
riodicity when the mutual statistical angle takes special 
values 9 = ±1/2, a Z 2 mutual CS theory must be en- 
forced on the lattice. Again, two ways of regularization 
may lead to different conclusions. In fact, the periodicity 
of the conventional CS theory under 9 — > 9 + 2 is also 
a very intricate issue. On the one hand, the perturba- 
tive RG expansion in terms of the statistical angle 9 in 
the continuum limit in Refs. |2i"|-|23j| does not have this 
periodicity. On the other hand, properly regularized on 
the lattice, the CS theory does have this periodicity [E8|. 
The two different regularization do lead to different con- 
clusions on the Quantum Hall transitions. It also leads 
to notorious Hall conductivity difficulty at v = 1/2 [p9| . 
Although we are unable to solve the critical behaviors 
of this combined Z 2 and U(l) gauge fields coupled to 
spinons and hc/2e vortices, we stressed the importance 
of putting the theory on the lattices to keep the periodic- 
ity of Z 2 mutual CS term and shed considerable lights on 
the physical picture and structure of the quantum tran- 
sitions driven by hc/2e vortices. We suggest that it may 
be possible to trace out the Z 2 gauge fluctuation on the 
lattice first to get an effective action without the discrete 
gauge fields, then take continuum limit of it and perform 
RG calculation. Further work is needed to sort out all the 



possible phases and to understand the critical behaviors 
of the phase transitions between these phases. 
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APPENDIX A: THE MOST GENERAL FORM 
OF THE GAUGE PROPAGATORS 

Adding the gauge fixing terms ^{{d^af) 2 + (<9 M a*) 2 ) 
to Eqj2^, we can find the gauge field propagators by in- 
verting the matrix 

\ n 2 (fc) e ^ A fc A n 3 (fc)fc(<W- h w L ) + b ^ 

(Al) 

The results are: 

(r > _ ak^ky n 3 1 k^k v 

ic*wv- fc4 + n 1 n 3 + n 2 ; fc l(V,y & ' 

= (G^)^ - Ji^-i) 

In the Landau gauge a = and putting k — 0, we 
recover the results calculated in Sec. III. 



APPENDIX B: APPLICATION TO STATIC 
DISORDERED VORTEX ARRAY 

In Ref. [0] , Ye pointed out that the long range Loga- 
rithmic interaction between vortices suppress the super- 
fluid velocity fluctuation, but does not affect the internal 
gauge field fluctuation. He concluded that the quasi- 
particle scattering from the random gauge field domi- 
nate over that from the superfluid velocity ( the Volovik 
effect). In this appendix, by using the formalism devel- 
oped in the main text, we provide additional evidence for 
this conclusion. 

For random pinned vortices, the static vortices do not 
feel the statistical gauge field a* in Eqjlj]. This corre- 
sponds to a mass term for this gauge field. Adding the 
mass term to the equation and integrating out a*, it gen- 
erates a Maxwell term for |30t| : 

C = ^(d,-ia^ a + ±(fp 2 

+ \(d» ~ *a^\ 2 + + \ft v - iA e Jfe^ x d„a x (Bl 
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In the above Eq. the averages over a*, 



a M , $ should 



be understood as quenched instead of annealed averages. 
As discussed in the main text, the long-range logarithmic 
interaction between vortices represented by the Maxwell 
term suppress the fluctuation of the superfluid velocity. 
The quasi-particle is moving in a long-range correlated 
random magnetic field and short-range correlated scaler 
potential, in consistent with the physical picture in Ref. 



APPENDIX C: HC/E VORTICES, SPIN-CHARGE 
SEPARATION AND MOTT INSULATOR 

In sec III, neglecting the charge fluctuation, we con- 
centrated on the effect of mutual C-S term. In this ap- 
pendix, we do the opposite, namely neglecting the mu- 
tual C-S term, but concentrate on the effect of charge 
fluctuation. This is realized by using neutral-like trans- 
formation, namely setting <pA = i>B — 0/2 which is 
single-valued in the resulting hc/e vortices. The double 
strength vortex and its stability was investigated in Ref. 
Jl2| , the associated spin-charge separation has been dis- 
cussed extensively in Refs. ]l3| , p^| . In this appendix, to 
contrast with hc/2e vortices discussed in the main text, 
we review their results, fill in some details and also get 
the new result on the expression of current operator 1^ 
in the superconducting phase. 

Putting (f>A = 4>b = 0/2, namely setting a M = in 
Eq.[fl] leads to: 



£ = VlA^r +V f p x T 3 + V2PyT 1 ]lpi a + (1->2,X 



+ f(^0' 



v) 

(Cl) 



c = i>hMa + - la^ 2 + v(\<t>\ 



(C4) 



Setting af = A e J f up to irrelevant term leads to 

£ = ^alMa + n tfd T (f> +\{d M - lA^ 2 + + • • 

(C5) 

where the second term is the Berry phase term and 
• • • means the irrelevant couplings between bosons and 
fermions. The above Eq. is essentially the same as 



Eq.Cl. Bosons and fermions are asymptotically decou- 
pled. The fermions (spinons) carry only spin 1/2, the 
bosons ( holons) carry only charge e. There is spin- 
charge separation. The long-range Coulomb interaction 
could be incorporated by adding the dynamic term for 
the time component of the gauge field \k\A {k)\ 2 [ f3l"| . 
Without the fermionic part, the action is the same as 
the superconductor-insulator transitions in dirty metals 
studied in [p^] , except here condensed is a charge e boson 
instead of a Cooper pair. 

In the underdoped regime, the vortex condenses < 
$ >= $o which generates a mass term for a M in Eq.C2 
5^o( a /-0t- Integrating out the massive a M leads to a 
Maxwell term for Afff = A^ — K~ 1 J fl from which we 
can identify the automatically conserved electric current 

i» = ^ o {dlJ»-d^d y j v ) @. 

In the superconductor phase, the vortex <!> is mas- 
sive, therefore can be integrated out, it leads to the 
old Maxwell term j^if^v) 2 - Integrating out a M gen- 
erates a mass term ^m^,{A e J^) 2 . From the mass term, 
we can identify the automatically conserved electric cur- 
rent 1^ = m$( J M )t. Note that although the quasi-particle 
electric current is not conserved itself, the total electric 



where the effective gauge field is A e J f = A M - K x J ll . current Z M = m$(5 jUI/ - ^r L )Jv{q) is conserved 
There is no internal gauge field a^, because double vor- 
tices do not scatter the quasi-particles. 

Performing duality transformation to vortex represen- 
tation leads to 



£ = ^7^a + |(^-^)$| 2 

1 



(C2) 



where the last ( Berry phase ) term can be absorbed into 
A e Jf by redefining Af/f — > A e Jf + ifiS^ . Comparing 
with Eq. 19, there is only charge fluctuation. There is no 
mutual C-S term, namely 9 = or any integer . 

The duality transformation to boson representation is 

C = 4>t<yM a + |(a M - m£)0| 2 + 1/(|0|) 
+ {{ft? + ia^Ma* - A** 16) + \f 2 u (C3) 

Integrating out the electric gauge field leads to 
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